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We have performed synchrotron radiation X-ray and neutron diffraction measurements on mag- 
netoelectric multiferroic CuFei_ ;z Ab02 (x — 0.0155), which has a proper helical magnetic structure 
with incommensurate propagation wave vector in the ferroelectric phase. The present measurements 
revealed that the ferroelectric phase is accompanied by lattice modulation with a wave number 2q, 
where q is the magnetic modulation wave number. We have calculated the Fourier spectrum of 
the spatial modulations in the local electric polarization using a microscopic model proposed by 
Arima [T. Arima, J. Phys. Soc. Jpn. 76, 073702 (2007)]. Comparing the experimental results with 
the calculation results, we found that the origin of the 2g-lattice modulation is not conventional 
magnetostriction but the variation in the metal-ligand hybridization between the magnetic Fe 3 + 
ions and ligand O 2- ions. Combining the present results with the results of a previous polarized 
neutron diffraction study [Nakajima et ai, Phys. Rev. B 77 052401 (2008)], wc conclude that 
the microscopic origin of the ferroelectricity in CuFei-^ AUO2 is the variation in the metal-ligand 
hybridization with spin-orbit coupling. 



PACS numbers: 75.80.+q, 75.25.+Z, 77.80.-e 

I. INTRODUCTION 

Magnetically induced ferroelectricity or electric con- 
trol of magnetic ordering has been intensively investi- 
gated since a colossal magneto-electric (ME) effect was 
found in TbMnOs^ Recent experimental studies have 
discovered a variety of ferroelectric magnetic compounds, 
which might be suitable for use in advanced ME de- 
vices. These materials are often termed ME multifer- 
roics. In order to develop practical device applications, 
it is essential to determine the microscopic mechanism of 
spin-polarization coupling in these systems. So far, there 
have been two different microscopic models for magnet- 
ically induced ferroelectricity. One is the spin-current 
model^ which predicts that two noncollinearly aligned 
neighboring spins Si and Si+i generate a local electric 
dipole moment p given by p oc e^t+i x (Si x Si+i), 
where e^j+i is the unit vector connecting the two spins. 
This formula predicts macroscopic uniform electric po- 
larization in a magnetic structure with cycloidal spin- 
components, and shows excellent agreement with the ex- 
perimentally determined relationships between magnetic 
structures and electric polarization in various transi- 
tion metal oxides, such as TbMnC>3j2 Tbi-^Dy^MnOa, 4 
Ni 3 V 2 8 , 5 MnW0 4 , 6 and CoCr 2 4 ^ The other micro- 



scopic model is the magnetostriction model, which pre- 
dicts ferroelectricity in collinear-commensurate magnetic 
structures. In this case, the local electric dipole mo- 
ment is given by C(r)(Si ■ Si+i), where C(r) is a con- 
stant dependent on the local crystal structure and the 
exchange interactions. The spontaneous electric polar- 
ization observed in the collinear-commensurate magnetic 
orderings in orthorhombic HoMnOs 8 may be attributed 
to this model. However, there are several ME multiferroic 
materials whose ferroelectricity cannot be explained by 
either of these models^ for example, delafossite multifer- 
roic CuFe02' 10111 ' 12 Hence, identification of microscopic 
spin-polarization coupling in these systems paves the way 
for a new design of multiferroic materials. 

The crystal structure of CuFeC>2 is shown in Figs. QJa)- 
(c). Owing to the geometrical frustration in the trian- 
gular lattice planes of the magnetic Fe 3+ ions, CuFe02 
exhibits various magnetically ordered phases i 13 ' 14 ! 15 The 
ground state of CuFe02 is a collinear-commensurate 4- 
sublattice (4SL) antiferromagnetic state. A spontaneous 
electric polarization emerging in the direction perpendic- 
ular to the hexagonal c axis was discovered in the first 
field-induced phase.™ Subsequent studie d 16 ' 17 revealed 
that the ferroelectric phase is stabilized even under zero 
field by substituting a small amount of nonmagnetic Al 3+ 
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FIG. 1: (Color online) (a) Crystal structure of CuFe02 with 
the hexagonal basis (qj,, bh, c; t ). (b)-(c) Relationship between 
the hexagonal basis and the monoclinic basis (a,b,c). (d) 
Schematic drawing of the proper helical magnetic structure 
in CuFei_ I Al a; 02. 



ions for magnetic Fe 3+ ions, as seen in Fig. [5J Re- 
cent neutron diffraction measurements in applied field 
revealed that the ferroelectric phase magnetic structure 
is an antiferromagnetically stacked proper helical mag- 
netic structure fii The wave vector is incommensurate. 
The helical axis is parallel to the hexagonal [110] direc- 
tion, as shown in Fig. [ljd). Hereafter, we refer to the 
ferroelectric phase as the ferroelectric incommensurate 
(FEIC) phase. Since neither the spin-current model nor 
the magnetostriction model is able to explain the ferro- 
electricity in proper helical magnetic ordering, we antic- 
ipate that another type of spin-polarization coupling is 
realized in CuFei_ :c AL r 02. 

Recently, Arima presented a theoretical consideration 
on the ferroelectricity in CuFei-^ALjC^ suggesting that 
the variation in the metal-ligand hybridization with spin- 
orbit coupling is relevant to the ferroelectricity in this 
system. Applying a microscopic theory derived by Jia 
et a/ . 25 ' 26 to a cluster model with a proper helical spin 
arrangement, Arima has predicted the following intrinsic 
features of the ferroelectricity in this system: 

Feature (i): the direction of the uniform polarization 
should be parallel to the helical axis of the proper helical 
magnetic ordering. 

Feature (ii): the spin helicity (i.e., right- or left-handed 
helical arrangement of spins) should correspond to the 



polarity of the uniform polarization. 

Feature (iii): There must be spatial modulations with 
wave numbers of 2q and 4g (where q is the magnetic 
modulation wave number) in the helical-axis components 
of the local electric polarization vectors. 

The results of recent polarized neutron diffraction mea- 
surements on CuFei_ K Al x 02 (x — 0.02) under an applied 
electric field show excellent agreement with features (i) 
and (ii)i^ However, the existence of spatial modulations 
of the local polarization (feature (hi)) has not yet been 
confirmed. Spatial modulations of the local polarization 
must result in lattice modulations, which can be observed 
by X-ray diffraction measurements. The existence of in- 
commensurate lattice modulations has been reported in 
the field-induced FEIC phase of undiluted CuFe02 by 
Terada et al^ and Ye et alM- However, application of 
a magnetic field readily induces slight distortion of the 
magnetic structure as well as the uniform magnetization 
components. Actually, higher harmonic magnetic reflec- 
tions have been detected in the field-induced FEIC phase 
of CuFe0 2 ^ As discussed by Terada et al. in Refl27l. 
this situation can easily cause additional lattice modu- 
lations through magnetostriction. In order to elucidate 
the existence of the lattice modulations predicted by the 
metal-ligand hybridization model (Feature (iii)), it is es- 
sential to remove the contribution of these field-induced 
lattice modulations. 

In the present study, we have performed synchrotron 
radiation X-ray measurements using a CuFei_ 2; AL I ;02 
(x = 0.0155) sample, which exhibits a FEIC phase un- 
der zero field. The present measurements reveal the ex- 
istence of 2g-lattice modulations in the zero-field FEIC 
phase. We have also investigated the magnetic-field de- 
pendence of the 2q lattice modulations by X-ray and neu- 
tron diffraction measurements under a field applied along 
the helical axis. In order to identify the origin of the 
observed 2q-lattice modulation, we calculate the Fourier 
spectrum of the spatial modulations in the helical-axis 
component of the local electric polarization using the 
microscopic model presented by ArimaJ^ This calcula- 
tion reveals that the lattice modulation does not originate 
from a conventional magnetostriction but rather from the 
variation in the metal-ligand hybridization between the 
magnetic Fe 3+ ions and the ligand O 2- ions. As a re- 
sult, we conclude that the microscopic mechanism of the 
ferroelectricity in CuFei_ 2 ,Al a; 02 is the variation in the 
metal-ligand hybridization with spin-orbit coupling. 



II. EXPERIMENTAL DETAILS 

A single crystal of CuFei_ a; Al :E 02 with x = 0.0155 was 
prepared by the floating zone technique^, and was cut 
into a disk shape. The Al-concentration was determined 
by chemical analysis. We performed X-ray diffraction 
measurements in zero field at the beamline BL46XU in 
SPring-8. The energies of the incident X-ray beams were 
tuned to 12 keV. The sample was mounted in a closed- 
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FIG. 2: (Color online) Schematic drawing of x-H-T phase 
diagram of CuFei-a; A\ x 02- 



cycle 4 He-refrigerator. The data obtained in this mea- 
surement are presented in sec. III. A. 

We have also performed X-ray diffraction measure- 
ments in zero field and in fields applied along the hexago- 
nal [110] (monoclinic b axis) at the beamline BL22XU in 
SPring-8. The sample was mounted in a horizontal- field 
cryomagnet whose maximum field is 6 T; the hexagonal 
(hhl) plane was used as the scattering plane. The data 
obtained in these measurements are presented in sec. Ill, 
B and C. 

Neutron diffraction measurements in a field applied 
along the hexagonal [110] axis were performed at the 
two-axis neutron diffractometer E4 installed at the 
Berlin Neutron Scattering Center in Hahn-Meitner In- 
stitute. The hexagonal (hhl) plane was selected as the 
scattering plane. Incident neutrons with wave numbers 
of 2.44 Awere obtained by a pyrolytic graphite (002) 
monochromator. An external magnetic field directed 
along the hexagonal [110] direction was provided by the 
horizontal field cryomagnet, HM-2, whose maximum field 
is 4 T. 

As described later, CuFei_ a; AL J 02 with x — 0.0155 
exhibits a symmetry-lowering structural transition from 
a rhombohedral structure to a structure with lower sym- 
metry (probably monoclinic symmetry, but possibly even 
lower symmetry). Taking into account the monoclinic 
lattice distortions found in the magnetically ordered 
phases of undiluted CuFe02 , 21 ' 22123 it is reasonable to 
employ a monoclinic basis in addition to the conventional 
hexagonal basis. The definitions of these bases are shown 
in Figs. [Tfa)-(c). We mainly employed the monoclinic 
notation. To distinguish between the two bases, the sub- 
script 'h' has been added to the hexagonal notation when 



referring to modulation wave numbers and reciprocal lat- 
tice indices. 



III. EXPERIMENTAL RESULTS 

A. Characterization of the crystal structures in the 
magnetically ordered phases 

Before discussing the lattice modulations in the ferro- 
electric phase, we present the symmetry-lowering struc- 
tural transition corresponding to the magnetic phase 
transitions in CuFei_ a; AL I .C)2 with x — 0.0155. At first, 
we should review the magnetic phase transitions in the 
x = 0.0155 sampled As shown in Fig. [2J this system 
exhibits three magnetically ordered phases in zero field. 
For all the phases, the magnetic modulation wave vec- 
tors are described by (0,q, |)[= (q h ,q h , §) fc ]. The high- 
est temperature phase and the intermediate phase are 
the oblique partially disordered (OPD) phase and the 
partially disordered (PD) phase, respectively. Both of 
them have collinear sinusoidal magnetic structures with 
incommensurate wave numbers. The wave number of 
the OPD phase, q OPB ~ 0.390, is independent of tem- 
perature, while that of the PD phase, g PD varies with 
temperature (0.404 < q PB < 0.430). The lowest temper- 
ature phase is the FEIC phase, whose modulation wave 
number is g FEIC = 0.414. 

Figures EHa)-(d) show the X-ray diffraction intensity 
maps around the reciprocal lattice position of (2,2,0)^ 
in the paramagnetic (PM) phase and the magnetically 
ordered phases. As seen in Fig. [3(d), a single peak as- 
signed as (2, 2, Q) h is observed at T = 19.2 K in the PM 
phase. In the PD and FEIC phases, this fundamental 
peak splits into several peaks, as shown in Figs. [3(a)- 
(b). This suggests that the threefold rotational symme- 
try along the Ch axis vanishes, resulting in a monoclinic 
(or even lower symmetry) lattice distortion in the PD 
phase and the FEIC phase. 

The monoclinic lattice distortions are also observed 
in the undiluted (x = 0.00) system , 21 > 22 in which the 
PD and FEIC phases show up as the thermally induced 
phase and the first field-induced phase, respectively. As 
discussed in the previous studies ; 21 ! 22 ' 23 these structural 
transitions are due to the bond order induced by the 
magnetostriction, which leads to lower-symmetry mag- 
netic orderings and lifts the macroscopic degeneracy of 
the magnetic states. 

On the other hand, splitting of the (2, 2, 0)/j peak was 
not observed in the OPD phase, as shown in Fig. [3(c), 
even though the magnetic structure of the OPD phase 
does not have threefold rotational symmetry along the Ch 
axis. 24 This suggests that the coherent bond order is not 
essential to the OPD magnetic ordering. Based on the 
fact that the OPD phase is never observed without non- 
magnetic substitution, it is reasonable to propose that 
the OPD magnetic ordering is stabilized by local symme- 
try breaking due to the site-random magnetic vacancies. 
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FIG. 3: The X-ray diffraction intensity maps obtained by the (H,K, 0)h reciprocal lattice scans around the reciprocal lattice 
position of (2, 2, 0) h at (a) T = 5.5 K, (b) T = 10.5 K, (c) T = 11.0 K and (d) T = 19.2 K. These intensity maps were measured 
in a warming process, (e) Temperature variations of the lattice constants a and b in the monoclinic notation, which were 
measured on heating, (f) Temperature variations of the intensities of (0, 3, 2) and (0, 2 + 2g PD , 2) reflections. The inset shows 
the diffraction profiles of the (0, 2 + 2q m , 2) reflections at typical temperatures. The data of (0, 2 + 2qPD, 2) are magnified 
by a factor of 20. 



Using the monoclinic basis, we can identify the split- 
ting peaks in the FEIC and PD phases, as shown in Figs. 
131(a) and[3fb), and we can also estimate the temperature 
variation of the lattice constants a and b in the x = 0.0155 
sample, as shown in Fig. [3Je). The differences between 
the lattice constants in the FEIC (or PD) phase and those 
in the PM phase in this system are comparable with those 
in the undiluted system i 21 ' 22 ' 23 These results suggest that 
the crystal structures of the FEIC and PD phases in the 
x = 0.0155 sample are almost the same as those in the 
undiluted system. 

In the present measurements, we observed incommen- 
surate superlattice reflection at the reciprocal lattice 
point of (0, 2.86, 2) in the PD phase. This reflection can 
be assigned as T even + (0, 2g PD ,0), where T even is the 
reciprocal lattice point of (H, K, L) with the condition 
'if + K = 2n' (n is an integer). We also observed com- 
mensurate superlattice reflection at the reciprocal lattice 
points of (0,3,2) in the FEIC phase. The temperature 
variations of the intensities of these reflections are shown 
in Fig. [3]^f). The incommensurate superlattice reflection 
in the PD phase and the commensurate superlattice re- 
flection in the FEIC phase are also observed in the undi- 
luted system. ) 22 i 23 i 27 and consistently explained by the 
magnetostriction model proposed by Terada et al£l 



B. Lattice modulations in the FEIC phase 

We now focus on the incommensurate lattice modu- 
lations in the FEIC phase. Figures Sfb) and 31(c) show 
the X-ray diffraction profiles of the (0, K, 0) and (0, K, A) 
reciprocal lattice scans at T = 4.5 K in the FEIC phase 
under zero field. As shown in Fig. 4(b) and its inset, a 



satellite reflection is found at the reciprocal lattice po- 
sition of (0, 2.172, 0). The intensity of this reflection is 
smaller than that of the (0, 4, 0) fundamental reflection by 
a factor 10~ 7 ~ 10~ 8 . This reflection is successfully iden- 
tified as being (0, 3 - 2g FEIC , 0). It should be noted that 
the experimental resolution in the present measurement 
can clearly distinguish the difference between the q FElc 
and q PD phases at T = 4.5 K (~ 0.43), although the high- 
temperature PD phase coexists with the FEIC phase un- 
der zero field cooling in the x — 0.0155 sample . 18 i 28 We 
confirmed the repetition of this reflection in reciprocal 
lattice space, as shown in the Figs. Sfd) and[U(e). Here- 
after, we refer to these reflections as t d m ± (0, 2g FEIC , 0) 
reflections, where T dd is the reciprocal lattice point of 
(H, K, L) with the condition L H + K = 2n + V. 

The T dd ± (0, 2g FEIC , 0) reflections are observed only 
in the FEIC phase. However, as shown in Fig. [U(f), 
the temperature variation of the integrated intensity of 
the r dd ± (0, 2<7 FEIC , 0) reflection is not proportional to 
that of the magnetic order parameter; specifically, the 
intensity of the (0, 3 — 2g FEIC , 0) reflection obviously in- 
creases around T = 6.5 K with increasing temperature, 
and rapidly decreases at the transition temperature from 
the FEIC phase to the PD phase, T = 8 K, while the 
FEIC magnetic order parameter (neutron diffraction in- 
tensity) monotonically decreases with increasing temper- 
ature, as seen in Fig. |U(f). 

No significant reflections were detected in the (0, K, |) 
reciprocal lattice scan, as shown in Fig. Etc). The 
satellite reflections assigned as (0, 2n ± g FE , A) and 
(0, (2n+ 1) ±g FEIC , i) were observed in the field-induced 
FEIC phase of undiluted CuFe0 2 i^ These satellite re- 
flections indicate the existence of a lattice modulation 
with a wave number lq FEIC in the field-induced FEIC 
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FIG. 4: (Color online) (a) Illustration of the reciprocal lattice (OKL) zone, [(b)-(c)] Diffraction profiles of the (b) (0K0) and 
(c) (0K^) reciprocal lattice scans in the FEIC phase. Contamination peaks from the Cu sample holder are masked by the 
shaded areas. The inset of (b) shows the magnification around the reciprocal lattice position of (0, 3 — 2g FEIC , 0). [(d)-(e)] The 



diffraction profiles of (d) (0, 3 + 2q ,0) and (e) (0, 5 — 2q ,0) reflections, (f) Temperature dependences of the integrated 
intensity of the (0, 3 — 2g FBIC , 0) reflection and the neutron diffraction intensity of the magnetic Bragg reflection corresponding 
to the FEIC magnetic ordering (taken from Ref. [la ), under zero field. 



phase. However, these reflections were not observed in 
the zero-field FEIC phase of CuFei-^ALjC^ with x = 
0.0155. This is consistent with the calculation of Ter- 
ada et al.^ which shows that the combination of the 
proper helical spin components and uniform magnetiza- 
tion components along the Ch axis induces lq-lattice mod- 
ulation through magnetostriction, and suggests that the 
lg-latticc modulation is not essential to the ferroelectric- 
ity in this system. 

The present results reveal that the 2g-latticc modula- 
tion corresponding to the satellite reflections at t^j ± 
(0, 2<7 FEIC , 0) exists in the zero-field FEIC phase. The mi- 
croscopic origin of the 2<7-lattice modulation is discussed 
in sec. IV. 



C. Hui, dependence of the 2g-lattice modulation 

We also surveyed the magnetic field dependence of 
the T dd i (0, 2<7 FEIC , 0) reflections by the X-ray diffrac- 
tion measurements under an applied field. As mentioned 
above, application of a magnetic field along the Ch axis 
readily leads to additional lattice modulations through 
magnetostriction. We thus applied a magnetic field along 
the helical axis of the proper helical magnetic structure, 
i.e. along the b axis. When a magnetic field is applied 
along this direction, the proper helical magnetic struc- 
ture deforms into a conical magnetic structure, as shown 
in Fig. EJd). In this configuration, the magnetic field 



(i/||h) does not modify the amplitude of the magnetic 
moment at each Fe 3+ site. Hence, the magnetostriction 
cannot induce additional lattice modulations. 

Figure[5Ja) shows the integrated intensity of the (0, 3 — 
2(j FEIC ,0) reflection at T = 4.5 K as a function of the 
magnetic field along the b axis (H\\ b ). The intensity of the 
reflection increases with increasing magnetic field. This 
result leads to the following two possibilities: that the 
amplitude of the 2g-lattice modulation itself increases 
with increasing magnetic field, or the fractions of the 
three magnetic domains, which reflect the threefold rota- 
tional symmetry of the original trigonal crystal structure, 
change with magnetic field intensity. 

We have also performed (0, K, 0) and (0, K, i) recipro- 
cal lattice scans at T = 4.5 K under H\\ b = 6 T. In both 
these scans, no additional (field-induced) reflections are 
detected. 



D. Neutron diffraction measurements under 
applied field along the b axis 

In order to elucidate the origin of the iJ||jj-variation of 
the intensity of the (0,3 — 2g FEI , 0) reflection, we per- 
formed neutron diffraction measurements under an ap- 
plied field along the b axis. We found that the magnetic 
Bragg reflections corresponding to the FEIC magnetic 
ordering increase with increasing magnetic fields Fig- 
ure EJb) shows the field dependence of the intensity of 
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FIG. 5: (Color online) (a)The _ff||i,-dependence of the inte- 
grated intensity of (0, 3 — 2g FEIC , 0) reflection observed in the 
X-ray diffraction measurements, (b) The ffy^-dependence of 
(I — 9h EIC i \ ~~ I™ 10 , \ )h magnetic Bragg reflection observed 
in the neutron diffraction measurements at T — 4.5 K. (d) The 
field dependence of the integrated intensity of (0, 3 — 2g FEIC , 0) 
reflection normalized by the data of the neutron diffraction 
measurements, (d) Schematic drawings of a conical structure 
with uniform magnetization component along the helical axis. 
[(e)-(f)] Schematic drawings of the fractions of three magnetic 
domains (e) under zero field, and (f ) under applied field along 
the b axis ([110];, axis). The directions of the filled arrows de- 
note the [001]h projections of the propagation wave vectors. 
The sizes of arrows qualitatively show the fractions of each 
domain. 



the (~ - g FEIC , 5 - if® 10 , §)/* magnetic Bragg reflection, 
where g FEIC [= ^g FEIC ] = 6.207. Since the change in the 
magnetic structure factor is considered to be negligible 
for magnetic fields of < < 4 T, 2 — this enhance- 
ment of the magnetic Bragg intensity indicates that the 
fraction of the FEIC magnetic domain increases with the 
propagation wave vector (g FEIC , <? FEIC , §)/,- This result 
supports the latter scenario for the enhancement of the 
intensity of the (0,3 — 2g FEIC , 0) reflection. A magnetic 
field along the b axis ([110]?, axis) should favor the mag- 
netic domain in which a more uniform magnetization 
component is induced along the direction of the mag- 
netic field. As a result, the fraction of the magnetic do- 
mains with a wave vector (g/,, <?/,, § )h is enhanced by the 
magnetic field along the [110]/, axis, and the fractions 
of the other magnetic domains with the wave vectors 
(-2g/,,g/,, §)/, an d {lh,—^1h,%)h should be reduced, as 
shown in Figs. [3£e) and[^f). The details of the magnetic 
domain distributions and the magnetic phase transitions 
in CuFei-^Al-rC^ (x — 0.0155) under applied field along 
the [110]/,, [110]?, and [001]/, directions will be presented 
in another paper i 2 ^ 

Figure^c) shows the intensity of the (0, 3 - 2g FEIC , 0) 
reflection normalized to the intensity of the magnetic 



Bragg reflection. In the magnetic field region of Hii^ < 
4 T, the amplitude of the 2g-lattice modulation does 
not exhibit remarkable change, and it decreases slightly 
above H\\ b = 3 T. 

IV. DISCUSSIONS 
A. Review of Arima's cluster model 

In this section, we discuss whether the microscopic ori- 
gin of the 2q-lattice modulation observed in the FEIC 
phase is the variation in the d-p hybridization or not. 
At first, we start from the Fe402 cluster model used in 
Ref. [H (see Fig. [B£a)). Applying the microscopic theory 
by Jia et all 25 ' 26 to the Fe-0 covalent bonding, Arima 
pointed out that the d-p hybridization between a mag- 
netic Fe 3+ ion and a ligand 2 ~ ion is slightly modified 
depending on the direction of the magnetic moment at 
the Fe 3+ siteJ^ As a result, the charge transfer from the 
magnetic Fe^ (i — 1,2,3,4) site to the neighboring Oj 
(J = 1, 2) site can be described by 

C* + AC cos 29., d , (1) 

where Co is the charge transfer in the PM phase, and 
AC is the constant dependant on the difference in the 
spin-orbit interaction between Fez and Oj, and the mag- 
nitude of the ordered magnetic moment at the Fe site, 
and so on. Qij is the angle between the direction of the 




FIG. 6: (Color online) (a) A Fe402 cluster used for the cal- 
culation of the local electric dipole moment in Reffl9l. (b) 
The magnetic structure in the FEIC phase with the a x b x 2c 
cell, (c) The phase differences in the magnetic modulations 
(A$™ 9 ) between the Fe 3+ sites. 
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Fej-Oj bond and the direction of the magnetic moment 
at the Fe^ site. The macroscopic uniform electric polar- 
ization is not induced by this term only. Arima argued 
that the many-body effect among the three covalent Fe- 
O bonds which are concerned with an 2 ~ ion, should 
be taken into account in Eq. (p}. As shown in Fig. EJa), 
an O 2- ion in this cluster (or CuFeC^) is surrounded by 
three neighboring Fe 3+ ions. Hence, the charge transfer 
between Fei and Oj should be slightly affected by the 
amount of charge transfer in the other two Fe-0 bonds 
which are concerned with the Oj site. For example, the 
covalency between Oi and Fei should be reduced when 
Fe2-Oi and Fe4-Oi bonds are more covalent. Taking ac- 
count of this many-body effect, the charge transfer from 
Fei to Oi can be described as follows^ 

C* [l + C"cos26 M ][l - aC (cos 20 2 ,i + cos20 4 ,i)], (2) 

where C = 1). a is the parameter representing 

the efficiency of the many-body effect, and is assumed to 
be small q« 1. Here, we apply to the cluster, a proper 
helical magnetic structure with a modulation wave num- 
ber q, whose helical axis is parallel to the Fei-Fe4 bond 
direction. The helical-axis component of the induced 
electric dipole moment at the Oi site, pb, should be pro- 
portional to the imbalance between the charge transfer 
from Fei to Oi and that from Fe 4 to Oi, 

Pb oc (1 + a)(cos 20i,i — cos 264,1) 

-aC"cos26 2 ,i(cos26i,i -cos20 4 ,i). (3) 

In the case of the proper helical magnetic structure, the 
value of cos 26;,j oscillates with a period of half the mag- 
netic modulation. As a result, the first term of the Eq. 
([3"|) gives the nonuniform polarization oscillating with a 
wave number of 2q (2g-modulation), and the second term 
of Eq. ([3]) gives the uniform polarization and the nonuni- 
form polarization oscillating with a wave number of Aq 
(Aq- modulation). It should be noted that one can de- 
rive Eq. (6) in Ref. [l9| from Eq.(j3j in this paper. The 
above formula of the local electric polarization suggests 
that the amplitude of the 4q-modulation is much smaller 
than that of the 2q- modulation. In the limit of a = 0, 
the uniform polarization and ^-modulation vanish. It is 
worth mentioning here that the same amount of electric 
dipole moment is induced at the O2 site, because of the 
symmetry of the proper helical magnetic structure. 



B. Calculation results in zero field 

We apply Eq. © to the magnetic structure in the 
FEIC phase, which has been determined by a previous 
magnetic structure analysis For simplicity, hereafter, 
we refer to the wave number and the wave vector in the 
FEIC phase as q and q[= (0, q, |)], respectively. In order 
to define the spin arrangement in the FEIC phase, we 
employ a x b x 2c cell, as shown in Fig. [BJb). This cell 



contains four Fe 3+ sites, and the fractional coordinates 
of these sites are: 

di = (0,0,0) 
d 2 = (0,1/2,0) 
d 3 = (0,0,1) 

d 4 = (0,1/2,1). (4) 

From the results of a previous magnetic structure anal- 
ysis, the spin components at di site are described, using 
the Cartesian coordinates shown in Fig. [DJc), as follows: 

Sf = fi x cos(27rq • (I + di) ~ fc) 
Sf = 

Sf = /i, «rin(27rg- (I + (5) 

fi x and fi z are the magnetic moments along the x and 
z axes, respectively. Since no significant ellipticity has 
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FIG. 7: (Color online) [(a)-(b)] 8 mas dependence of the dif- 
ference between the phases of (a) 2q- and (b) 4g-modulations. 
[(c)-(d)] The phase differences of the (c) 2q- and (d) Aq- 
components in the spatial modulations in the local electric 
polarization. 



been detected in the magnetic structure analysis for zero 
field (and for relatively low fields, -f^n [001] ^ < 4T) on the 
x = 0.0155 sample,— we assume at this stage that the 
magnetic structure has no ellipticity, i.e., fi x = fi z = fi. 
4>i is the relative phase shift at the di site; specifically, 



H - 02 | = 103 ~ 04 1 = S mag 
!>1-0 3 | = |02 - 4 |=O 



(6) 



Although S mag was determined to be ~ 76° (~ nq) by 
magnetic structure analysis^ we treat 8 mag as a param- 
eter in the following calculation. Here, we define the 
phase difference of the magnetic modulations between 
the di site and the dj site as A$™ aff . Using this nota- 
tion, the relationship among the phases of the magnetic 
modulations on the four Fe 3+ sites is summarized in Fig. 
He). 

We numerically calculated the spatial modulations of 
Pb at each oxygen site using Eq. (J3J) . The system size is 
set to be a x Nb x 2c. Since the magnetic propagation 
wave number along the b axis is incommensurate, N is set 
to be a large number (typically, N ~ 100). We also calcu- 
lated Fourier spectra for their nonuniform components. 
Here, we define the phase difference in the rig-modulation 
(n = 2, 4) between the oxygen site neighboring di site and 
that neighboring dj site, as A$™?, in the same manner 
as that of the magnetic modulations. The results of the 
calculations for A<3?^- and A<&$ are shown in Figs. Ha) 
and [7(d) . For the 2q- modulations, there is a finite phase 
shift, 25 mag , between the oxygen site neighboring the d\ 
(d^) site and that neighboring the di (d^) site. For the 
4q-modulation, there are no phase shifts along the a axis. 
For both of the 2q- and 4g-modulations, there are also no 
phase shifts along the c axis. These results are summa- 
rized in Figs. [7(c) and [7(d). 

The intensities of the calculated Fourier spectra for 
the spatial modulations of pb, \P(Q)\ 2 , where Q is a 
vector in the reciprocal lattice space, are mapped onto 
the reciprocal lattice space, as shown in Figs. [8(a) 
and [5(b). Since the d-p hybridization model predicts 
that the local electric polarization arises from the lo- 
cal imbalance of the Fe-0 bond-covalency, it is reason- 
able to consider that the spatial modulations of the lo- 
cal electric polarization result in spatial modulations of 
the local atomic displacements of the oxygen ions. We 
thus assume that the local displacement of the oxygen 
ion is proportional to the magnitude of the local po- 
larization vector. In this case, \P(Q)\ is proportional 
to the structure factor in the X-ray diffraction mea- 
surements. In the case of 6 mag = 0, tt, satellite peaks 
corresponding to the 2q- and 4q-modulations appear at 
T even ± (0, 2q, 0) and r even ± (0, 4g, 0), respectively. In the 
case of 5 mag ^ 0, w, satellite peaks corresponding to the 
2(j-modulation appear, in addition to the above peaks, at 

76°, the ratio of 



T ol id ± (0, 2q, 0). By substituting 5. 

the intensities of the spectra, \P(T even ± (0, 2q, 0))| 2 and 
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FIG. 8: (Color online) Fourier spectrum mappings of the spa- 
tial modulations in the local polarization, which are calculated 



using (a) the d-p hybridization model with 5„ 



0,7T, (b) 



that with 8m a g = 76°, and (c) the magnetostriction model 
with finite ellipticity. The size of a symbol, except for T even 
and T dd, qualitatively shows the intensity of the spectrum 
|P(Q)| 2 , which is normalized by the intensity of the most 
strongest peak in the spectrum. The black arrows denote the 
position of the satellite reflections observed in the zero-field 
FEIC phase. 



\P(r dd ± (0, 2q, 0))| 2 is calculated to be 
\P(t 

even 

±(0,2g,0))| 2 



|P(r ocW ±(0,2 g ,0))| 2 



= 0.062. 



(7) 



This suggests that in X-ray measurements, the satellite 
reflections corresponding to the 2(7-modulation should 
be mainly observed at the reciprocal lattice position of 
r d(i± (0, 2q, 0). This shows good agreement with the re- 
sults of the present measurements. Although reflections 
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at r even ± (0, 2q, 0) and T even ± (0,4q, 0) were not ob- 
served in the present measurements, this can be ascribed 
to the S /N of the present measurements, because the in- 
tensities of those reflections are estimated to be much 
smaller than that of r Q dd ± (0, 2q, 0) reflections. 

It should be noted that this result does not change 
in the parameter region of < a < 1, although we 
assumed a = 0.1 in the above calculation. This is 
because the value of a affects only the ratio between 
the amplitude of the 2g-modulation and that of the 
Aq- modulation (and uniform polarization); specifically, 
\P(r even ± (0, 4q, 0))\ 2 /\P(T odd ± (0, 2q, 0))| 2 is in the or- 
der of ~a 2 C" 2 (< 1) 

It is worth mentioning here that neither the magne- 
tostriction model nor the spin-current model explains 
the satellite reflections at T dd ± (0, 2g,0), even if the 
magnetic structure has finite ellipticity. The spin-current 
term, e^j+i x (Si x S^ + i), does not produce the helical- 
axis component of the local polarization vector, in the 
proper (or elliptic) helical magnetic ordering in this sys- 
tem. Therefore, the lattice modulations induced by this 
term do not contribute to the satellite reflection in the 
reciprocal lattice position of (0, K, 0). The magnetostric- 
tion term, (Si ■ Si+i), produces the spatial modulation 
with a wave number of 2q in the helical-axis compo- 
nents of the local polarization vectors, when the mag- 
netic structure has finite ellipticity— However, as dis- 
cussed by Terada et ai.j21 the satellite reflections corre- 
sponding to the magnetostriction-induced 2q-modulation 
are observed only at the reciprocal lattice position of 
T even ± (0, 2q, 0), regardless of the value of 5 mag (see Fig. 
IHfc)). Therefore, we conclude that the origin of the 2q- 
lattice modulation observed in the zero-field FEIC phase 
is the variation in the metal-ligand hybridization with 
spin-orbit coupling, which corresponds to the feature (iii) 
mentioned in the introduction. 

We should mention the non-monotonic temperature 
variation of the intensity of the r a dd ± (0, 2q, 0) reflec- 
tion observed in the present measurement (see Fig. HIT)). 
The d-p hybridization between the Fe 3+ ions and the 
2 ~ ions must depend on the length of the Fe-0 bonds 
as well as the angle between the Fe-0 bond direction 
and the magnetic moment at the Fe site. We therefore 
anticipate that the temperature variations of the lattice 
constants and the fractional coordinates of the oxygen 
sites affect the temperature variation of the intensity of 
the T dd ± (0, 2q, 0) reflection. The precise determina- 
tion of the structural parameters by high-resolution X- 
ray and/or neutron diffraction measurements are desir- 
able for further clarification. 



C. Calculation results for an applied field 

When a magnetic field is applied along the b axis, 
the magnetic structure is expected to be a conical struc- 
ture. We also numerically calculated the Fourier spectra 
for the nonuniform polarizations induced by the coni- 
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FIG. 9: (Color online) (a) 8 dependence of the intensities of 
the Fourier spectra for the spatial modulations in the local 
polarization. To enhance the visibility, the data of 3q and 
4q are magnified by a factor of 200. (b) Fourier spectrum 
mapping of the spatial modulations in the local polarization, 
which are calculated with the d-p hybridization model with 
Smag = 76°, a — 0.1 and 9 = 10°. The size of the symbols, 
except for T even and T dd, corresponds to the intensity of the 
spectrum |P(Q)| 2 , which are normalized using the intensity 
of the most intense peak in the spectrum. The black arrows 
denote the positions of the satellite reflections observed in the 
present measurements. 



cal magnetic structure. Figures [Ha) and E{b) show the 
calculation results with the parameters of a = 0.1 and 
Smag = 76°. 9 is the angle defined in the Fig. [5ja). With 
increasing 9, the intensities of the Fourier spectra for 2q- 
and 4g-modulations decrease. In the region of < 6 < 
lq- and 3q-modulations emerge from the first and second 
terms of Eq. ([3]), respectively. In reciprocal lattice space, 
the satellite peaks corresponding to the I ^-modulations 
appear at r even ± (0, q, |) and t 0(M ± (0, q, |), and the 
satellite peaks corresponding to the 3g-modulations ap- 
pear at T even ± (0, 3q, i;) and T a dd ± (0, 3q, ^), as shown 
in Fig. Mb). 

According to the results of recent magnetization 
measurements^ the canting angle 9 at Hn = 4T is 
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roughly estimated to be 5° ~ 10°. In the region of 
6 < 10°, \P(r odd ± (0, 2q, 0))| 2 decreases slightly with in- 
creasing 8, and the intensities of the other Fourier spectra 
are quite small compared with \P(T dd ± (0, 2q, 0))| 2 , as 
shown in Fig. [9ja). Hence, it is expected that the inten- 
sity of the T dd ± (0,2g,0) satellite reflection decreases 
slightly with increasing applied magnetic field along the 
b axis. This is consistent with the results of the present 
measurements. Although we could not survey the region 
of 9 > 10°, where remarkable changes in the amplitudes 
of the lattice modulations are expected, these results im- 
ply that the d-p hybridization model still works for a 
finite magnetic field along the b axis. 

V. SUMMARY 

In summary, we have performed synchrotron radia- 
tion X-ray and neutron diffraction measurements under 
zero field and applied field, on the delafossite multiferroic 
CuFei_ 2 ,Al a; 02 with x = 0.0155, in which the ferroelec- 
tric phase shows up under zero field. 

We found that the threefold rotational symmetry along 
Ch axis vanishes in the PD and FEIC phases. Although 
the present result does not determine the symmetry of 
the crystal structure in the PD and FEIC phases, it is 
reasonable to propose that a monoclinic lattice distor- 
tion occurs in the FEIC and PD phases in CuFei_ 2 ,Al a; 02 
with x = 0.0155, because the monoclinic lattice distor- 
tion is observed in the PD and FEIC phases in undiluted 
CuFe02 J 21l22 ' 2:3 On the other hand, monoclinic lattice 
distortion was not observed in the OPD phase, which 
never shows up without nonmagnetic substitution. This 
implies that the OPD magnetic ordering is stabilized not 
by the coherent bond order due to magnetostriction, but 
by the local symmetry breaking due to site-random mag- 
netic vacancies. 



In the FEIC phase, we found satellite reflections iden- 
tified as T dd ± (0, 2g FEIC , 0) under zero field. This indi- 
cates that the FEIC phase is essentially accompanied by 
2<7-lattice modulations. We have calculated the Fourier 
spectra of the spatial modulations in the local electric po- 
larization using the experimentally determined magnetic 
structure in the FEIC phase and the microscopic model 
presented by ArimaJ^ Comparing the experimental re- 
sults with the calculated Fourier spectrum revealed that 
the origin of the 2<7-latticc modulation is not conventional 
magnetostriction but the variation in the d-p hybridiza- 
tion between the magnetic Fe 3+ ions and the ligand 2 ~ 
ions. Combining the present results with the results of 
a previous polarized neutron diffraction study^ we con- 
clude that the microscopic origin of the ferroelectricity in 
this system is the variation in the d-p hybridization with 
spin-orbit coupling. 
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